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The string susceptibility exponents of dynamically triangulated two dimensional surfaces with sphere and torus 
topology were calculated using the grand-canonical Monte Carlo method. We also simulated the model coupled 
to d-Ising spins (d=l,2,3,5). 



1. Introduction 

For 2-dimensional quantum gravity the parti- 
tion function with a fixed topology is approxi- 
mated by 



Z(A) « A 1 ~ i e 
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for large A limit, where A is the 2-dimensional 
volume (area), A c is the cosmological constant, 
and 7 is the string susceptibility exponent. The 
exponent 7 is known to be 
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by analytical theories [|l]|2j] , where h is the number 
of handles of the 2-dimensional surface, and c is 
the central charge of matter fields coupled to the 
surface. When the dynamical triangulation (DT) 
recipe Q is adopted to generate random surfaces, 
it becomes possible to simulate the discretized 2- 
dimensional gravity numerically. 

In order to measure the exponent 7 with 
any topology, we introduce the grand-canonical 
Monte Carlo method j|. In this method the fi- 
nite size effect || enters through the total area 
and becomes controllable in reasonably large size 
simulations. Other method such as the MINBU 
(minimum neck baby universe) algorithm Q have 
been successful in the measurement of 7 of sphere 
topology for c < 1. However, in the MINBU 



algorithm measuring 7 for the case of a higher 
genus is more complicated than for a sphere, be- 
cause the minimum neck baby universe must not 
be topologically non-trivial and the baby universe 
distribution describing two topologically different 
7 must be separated independently |Q. Eq.(||) 
shows that 7 becomes complex for 1 < c < 25. 
However, in a numerical simulation by dynamical 
triangulation the partition function is kept real 
for any value of c. Then we evaluate the expo- 
nent 7 for c > 1 using the grand-canonical Monte 
Carlo method. 

2. Model and Method 

2.1. Model 

We start from the discretized Euclidean model. 
The discretized partition function with a fixed 
topology is given by 



Z(X) = J2 W{T,N 2 )e- XN *Z a , 
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where 
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Here W(T, N2) is the symmetry factor, T refers 
to a specific triangulation, Ni is the total num- 
ber of i-simplices, A is a parameter, and /3 is a 



coupling strength of Ising spins. We adopt the 
model excluded "tadpole" and "self-energy" dia- 
grams. There are d kinds of Ising spins on a 
vertex, which carry c = i central charge in total. 

2.2. Measurement of 7 

For the large- N2 limit the partition function is 
written as 
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Then, the distribution P(N 2 ) is given by P(N 2 ) oc 

N]- 3 e~ (x ~ K)N2 . From 

lnP(iY2) -(7-3)-(A-A c )-^-+ CO n^ (6) 
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if we tune A close to A c within the precision | A — 
A c I <C ^j^r > we can evaluate the exponent 7 by 
plotting \nP(N 2 ) as a function of lniV2 within 
the region iV 2 mm < iV 2 < N 7 2 nax . 

2.3. Grand-canonical simulation 

We explain the method used to generate DT 
surfaces in the case of pure gravity (c = 0). Sup- 
pose that we wish to increase the number of real 
vertices from N to N + 1; to do this we choose 
two dual links belonging to the same dual loop L it 
and insert a "propagator" dual link joining the 
two dual links, thus creating an additional dual 
loop Ln +i- For the inverse move, No — ► No — 1, 
we choose a dual link at random, and remove it, 
resulting in two dual loops Li separated by the 
dual link merging into one. The numbers of pos- 
sible increase and decrease moves from configura- 
tion A are 
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where {q{\ is the set of coordination numbers of 
a configuration A. 

In these moves the ergodic property is automat- 
ically satisfied. However, we must treat the de- 
tailed valance carefully || , since the total number 
of possible moves depends on N$ and The 
detailed valance equation, that the move should 
satisfy, between configurations A and B, is given 
by 
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where Pa->b is the transition probability from 
configurations A to B, Wa is the Boltzmann fac- 
tor determined from eq.(||), and tia is the total 
number of possible moves from configuration A. 



For the case c 



we can easily extend the 



method by modifying the total number of moves 
for including the freedom of spins as 



n A 



N 

E 



2 2 



(9) 



For the spin variables, Monte Carlo trials with 
Wolff's cluster algorithm [gj are performed after 
some geometrical moves. 

3. Results and Discussions 

We measured the string susceptibility exponent 
7 of the 2-dimensional surface (h = 0, 1) coupled 
to d-Ising spins (d = 0, 1, 2, 3, 5) using the grand- 
canonical Monte Carlo method. When Ising spins 
are put on the surface, we determine the critical 
temperature f3 c from the peak of magnetic sus- 
ceptibility for each combination of h and c. Es- 
pecially in the case of a torus (h = 1) the experi- 
mental values P c were quite close to the theoreti- 
cal value (P c w 0.2163). One sweep is defined by 
2000 x 2 d trial moves and 5 cluster operations of 
Ising spin. 

In order to obtain P(N 2 ), we measure the num- 
ber of triangles N 2 10 5 times every 5 sweeps 
within the region 1000 < N 2 < 3000. Our results 
are summarized in Table The errors indicated 
in parentheses for the last two digits are only sta- 
tistical errors estimated from the least squares 
fits. The systematic errors are estimated to be 
around 0.03 - 0.02. 

For the pure gravity (c = 0), Table |l| shows 
that both the exponent 7 and the cosmological 
constant A c are good agreements to the theoret- 
ical values (7 = -0.5, A c = |lnf^ « 1.12467). 
Furthermore, in Figs J^,^ the exponents 7's are re- 
produced closely to the theoretical predictions as 
long as c < 1. This proves our expect that the 
grand-canonical method is almost free from the 
finite size effects. Since we could check that our 
method is able to treat torus topology, it is easy 
to expand to the case of a higher genus. 



Table 1 

String susceptibility exponent 7 of sphere and torus 



Sphere (h = 0) Torus (h = 1) 





7 


A c 




7 


A c 




c = 0.0 


-0.517(27) 


1.12467 




2.014(24) 


1.12469 




c = 0.5 


-0.395(27) 


1.53716 


0.227 


1.995(23) 


1.52890 


0.2163 


c= 1.0 


-0.009(26) 


1.94848 


0.226 


1.971(22) 


1.93335 


0.2163 


c= 1.5 


0.112(27) 


2.35843 


0.225 


1.980(23) 


2.33841 


0.2163 


c = 2.5 


0.238(27) 


3.17513 


0.223 


2.018(25) 


3.14950 


0.2163 
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Figure 1. String susceptibility exponents 7's of a 
sphere with various central charge c's. Solid line 
represents eq.(2) with h = 0. 
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